Ray-based description of normal mode amplitudes in a range-dependent waveguide 
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An analogue of the geometrical optics for description of the modal structure of a wave field in 
a range-dependent waveguide is considered. In the scope of this approach the mode amplitude is 
expressed through solutions of the ray equations. This analytical description accounts for mode 
coupling and remains valid in a nonadiabatic environment. It has been used to investigate the 
applicability condition of the adiabatic approximation. An applicability criterion is formulated as a 
restriction on variations of the action variable of the ray. 
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I. INTRODUCTION 

We consider a method for analysis of the modal struc- 
ture of a scalar wave field in a two dimensional range- 
dependent waveguide. This method derived in Refs. [1- 
3] (see also [4]) is based on projecting the ray representa- 
tion of the wave field onto normal modes and evaluating 
the corresponding integrals using the stationary phase 
technique. As a result the mode amplitude is expressed 
through parameters of ray trajectories. This approach es- 
tablishes a relationship between the ray and mode repre- 
sentations of the wave field in a range-dependent waveg- 
uide and provides an approximate solution to the mode 
coupling equation. 

In describing ray trajectories we apply the Hamilto- 
nian formalism taken in terms of the action-angle vari- 
ables [5, G]. It is shown that each mode is formed by 
contributions from rays whose action variables up to a 
multiplicative constant are equal to the mode number. 
In the nonadiabatic environment the action varies along 
the ray path. Our objective in this work is to find out 
how the nonadiabaticity of rays manifests itself in varia- 
tions of mode amplitudes. In the nonadiabatic waveguide 
an initially excited single mode breaks up into a group of 
modes. We argue that the width of the group is propor- 
tional to the spread of ray actions at the given range. 

The formula for the mode amplitude, valid in the nona- 
diabatic waveguide, has been used to formulate an appli- 
cability condition of the adiabatic approximation. In con- 
trast to traditional approaches for solving this problem 
[7, 8] we proceed not from the mode coupling equation 
but from its solution (albeit approximate). An applica- 
bility criterion discussed in the present paper takes into 
account accumulation of errors with range. 

All results are obtained for the wave field that obeys 
the parabolic equation. This assumption is not necessary: 
the results are readily generalized to the field governed 
by the Helmholtz equation. We use the parabolic equa- 
tion approximation for two reasons. First, unlike the 
Helmholtz equation it can be easily solved numerically 
which is very important for testing our approximate for- 
mulas. Besides, the parabolic equation formally coincides 
with the Schrodinger equation and therefore formulas de- 
rived in this work may be applied in quantum mechanics 
to study the motion of a particle in a potential well with 



time-dependent parameters. 

In the present paper we investigate general proper- 
ties of wave propagation in a range-dependent waveguide 
by the example of the underwater sound channel. This 
choice is caused by interests of the authors and by the 
fact that the topics addressed here are of importance in 
underwater acoustics [7-10]. 

The paper is organized as follows. Section II provides 
a brief description of the ray representation of the wave 
field. Main relations of the Hamiltonian formalism ex- 
pressed in terms of both position-momentum and action- 
angles variables are presented. In Sec. Ill we consider 
the mode representation (in the WKB approximation) of 
the wave field and the mode coupling equation. We also 
outline a traditional method for estimating the validity 
region of the adiabatic approximation by using the per- 
turbation theory for the mode coupling equation. In Sec. 
IV formulas expressing mode amplitudes through param- 
eters of ray trajectories are derived. Section V focuses on 
connection between variations of the action variable and 
mode coupling. It is shown how a standard ray trac- 
ing can be used to estimate the number of interacting 
modes. An applicability criterion for the adiabatic ap- 
proximation is formulated as restriction on the spread of 
ray actions. These results are verified and illustrated us- 
ing numerical examples. The conclusions are summarized 
in the final section. 



II. RAY REPRESENTATION OF THE WAVE 
FIELD 

A. Parabolic equation approximation 

Consider a monochromatic wave field at a carrier fre- 
quency / in a two-dimensional acoustic waveguide with 
the sound speed c being a function of depth, z, and range, 
r. We shall assume that at each range point the sound 
speed profile has a single minimum. The wave field u(r, z) 
is governed by the Helmholtz equation [7, 8, 11] 



1 d du d 2 u 

dz 2 c 2 (r,z) 



n 2 
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(1) 



where f2 = 2nf. Select a reference sound speed Co such 
that |c(r, z) — cq\ << cq (in underwater acoustics such 



2 



choice of Co is always possible) and assume that graz- 
ing angles of the propagating waves are small. Then it 
is convenient to introduce an envelope function v(r, z) 
connected to u(r, z) by an expression 



ikr 



(2) 



with k = 0/co- This function is approximately described 
by the standard parabolic equation [7, 11, 12] 



2ik£ + ^ - 2#Uv = 0, 
or dz z 



where 



(3) 



the wave field. In the case of a point source the wave field 
is determined by Eq. (3) with an initial condition 



v(0,z) = S(z - z ). 



(10) 



Then all rays escape the same point (0, zq) and each tra- 
jectory can be labelled by the ray starting momentum p . 
The amplitude 



A = 



k 



2iri \dz/dpo\ 



(11) 



where /j, is the Maslov index or integral number of times 
that the ray passes through caustics [14]. 



U(r, z) = 



1 - 



c 2 (r, z) 



(4) 



Notice that Eq. (3) formally coincides with the time- 
dependent Schrodinger equation. In this analogy r, k^ 1 , 
and U(r, z) play roles of time, Planck's constant, and 
potential, respectively. 



C. Action-angle variables 

Since ray trajectories in the waveguide are oscillating 
curves, their analysis may be simplified by using the so- 
called action-angle variables [5, 6]. In mechanics the lat- 
ter are often applied to study oscillations of particles in 
potential wells. 



B. Hamiltonian formalism in terms of 
momentum-position variables 

In the geometrical optics approximation the wave field 
is formed by contribution from all rays (eigenrays) ar- 
riving at the observation point. The ray trajectory is 
defined by the Hamilton equations [6, 12] 



dz dH dp 
dr dp ' dr 

with the Hamiltonian 

H = p 2 /2 + U. 



dH 

dz 



(5) 



(6) 



The momentum p is connected to the ray grazing angle 
X through the relation 



p = dz/dr = tan%. 



(7) 



A contribution to the wave field from an eigcnray is 
presented in the form 



Ae 



IkS 



(8) 



where S and A arc the ray eikonal and amplitude, respec- 
tively. The eikonal S is an analog to Hamilton's principal 
function in classical mechanics and it is given by an in- 
tegral 



/ 



(pdz — Hdr) 



(9) 



running over the ray path [5, 12, 13]. An explicit expres- 
sion for the ray amplitude depends on the source exciting 



1. Range-independent waveguide 

First, consider a range- independent waveguide where 
the Hamiltonian H remains constant along the ray path 
(analog to the energy conservation law in mechanics). 
The action variable I is defined by the integral [5, 6] 



/= — 

2tt 



pdz = — 

7T 



dz^2 [H-U(z)], (12) 



where z mm and z max are the lower and upper ray turning 
depths, respectively, satisfying the condition U(z) = H. 
Equation (12) defines the "energy" H as a function of 
the action /. The ray trajectory is a periodic curve whose 
period along the r-axis (the cycle length), D, and the an- 
gular frequency of spatial oscillations, u), are determined 
by the relation 



dH _ _ 2tt 
~dl ~D' 



(13) 



Take one cycle of the unperturbed ray path that begins 
at the minimum of the trajectory. At the first half-cycle 
the canonical transformation, 



p = p(i,6), z = z(i,e), 



(14) 



from the position- momentum, (p, z), to the action-angle, 
(1,0), variables is determined by the equations [5, 6, 15] 



P = 



dG 

~dz~' 



dG 



with 



G{I,z)= ( dz^/2 [H(I) - U(z)} 

J 2 m in 



(15) 



(16) 
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being a generating function of the canonical transforma- 
tion. At this half-cycle the angle variable 9 varies from 
to 7r. The transformation is continued to the next half- 
cycle (tt < 9 < 2n) using relations p(I, 9) = —p(I, 2ir — 9), 
and z(I, 9) = z(I, 2tt — 9). The ray equations in the new 
variables take the trivial form 



dJ_ dH 



d0 dH ,n 



Note, that the so defined angle variable 9 varies from 
to 27r at a part of the trajectory beginning at one min- 
imum and ending at the next one. To make the angle 
variable continuous, its value should be increased by 2n 
at the beginning of each new cycle. Both functions in Eq. 
(14) are periodic in 9 with period 2tt. Correspondingly, 
an arbitrary function F(z) expressed through the action- 
angle variables becomes a periodic function of 9 and can 
be expanded in a Fourier series 



F(z) = ^F v (I)cosv6. 



(18) 



v=0 



The coefficients F v (I) are expressed analytically only for 
a few special examples of U{z) and F{z). But numerical 
evaluation of these coefficients can be easily performed 
using a standard ray code. If we denote the ray trajectory 
with the action variable I by z(r), then 



D(I) 



fD(I) 

x dr F(z(r min + r)) cos (w(/)(r min + r)) , (19) 

Jo 

where go = 1, Qu = 2 for v > 0, and r m ; n is a range 
corresponding to a minimum of the ray path. 



2. Range-dependent waveguide 

In a range-dependent waveguide canonical transfor- 
mations from (p, z) to (7, 9) variables arc different at 
different range points. At each particular range r the 
transformation is defined using the so-called reference 
waveguide: an imaginary range-independent waveg- 
uide whose cross-section coincides with that of a real one 
at the given range. The connection between pairs (p, z) 
and (7, 9) at range r is determined by formulas present 
in Sec. II C 2 that should be applied in a corresponding 
reference waveguide. Since the functions p(I,6), z(I,0), 
G(I,z), H(I), co(I), and F V (I) become different at dif- 
ferent ranges they acquire an additional argument r that 
labels the reference waveguide. 

The canonical transformation is determined by the 
equation [5] 



dS = pdz - Hdr = dG - OdI - Hdr, 



(20) 



where H is a new Hamiltonian. According to this relation 
H(I,e,r) = H(I,r)+A(I,0,r), (21) 

with 

dG(I,z,r) 



A(I,9,r) = 



dr 



(22) 



z=z(I,8,r) 

The Hamilton equations retains their canonical form 
dl _dB__ _dA dB__dH__ dA 



where 



w(I,r) = 



dH{I,r) 
dl ' 



(24) 



The ray cikonal in the new variables can be found by 
integrating Eq. (20). After some algebra we get 



S= f {IdO - (H + A) dr) 
Jo 



-G(z(0),I(0),0) + G(z(r),I(r),r) 



-9(0)1(0) - 9(r)I(r). 



(25) 



Initial and final values of the angle variable, 9(0) and 
9(r), present in the right hand side should be taken mod- 
ulo 2ir. 

The term A can be expressed through coefficients of an 
expansion 



dU(r,z) 



dr 



= ^ V v (I,r) cos 6 (26) 

z=z(I,6,r) „ = o 



analogous to Eq. (18). First, notice that 

(dH(I,r) dU(r,z) 



dr 



dr 



dz 



y / 2[H(T^)~ THr^)]' 



(27) 



where ± denotes the sign of p at a current half-cycle. 
From Eq. (12) it follows that 

8H(I,r) = q,(J, r) 
dr 7r 



dU(r,z) 



dz 



dr ^2 [H(I,r)-U(r, z) 



(28) 
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The right hand side, formally, can be considered as an 
integral over a half-cycle of the ray path in the reference 
waveguide. Then 



;(/, r)dz 



and 



y/2 [H(I,r)-U(r, zj] 
dH(I,r) 



dO 



dr 



= V (I,r). 



(29) 



Similarly it can be shown that 

dU(r, z) dz 



I 



dr y/2[H(I,r)-U(r,z)] 



Combining Eqs. (27), (29), and (30) we arrive at 



^ oo 

A = --E^) 



sin (vff) 



(31) 



Note that the coefficients V V (I, r) can be calculated using 
Eq. (19) in a reference waveguide corresponding to range 
r with F replaced by V(r, z) — dU(r, z)/dr. When evalu- 
ating the right hand side of Eq. (19) one should integrate 
over the ray path computed in the reference waveguide 
but not in the real one. In this calculation the argument 
r of function V(r, z) must be considered as a constant 
labelling the reference waveguide. 

The adiabatic approximation for rays is valid if the 
sound speed c is a so slow function of range that A be- 
comes negligible. Then the action variable I remains 
constant along the ray path [5, 7]. 



III. MODE REPRESENTATION IN THE WKB 
APPROXIMATION 

A. Eigenfunctions and eigenvalues 



4> s = 7r and (f>b — 0. In case when the upper (lower) mode 
turning points lie within a water bulk we have <p s = —ir/2 
(4>b = — 7r/2). The eigenfunctions are orthogonal and 
normalized in such a way that 



J dz ipmtpn = 5„ 



(33) 



In the WKB approximation the eigenvalue of the m-th 
mode is 

H m = H(I m ), (34) 
where I m satisfies the quantization rule [7, 16] 



kl m = to - 



2tt 



(35) 



The quantity I m can be treated as an action variable 
associated with the m-th mode: if only this mode is ex- 
cited then all rays have the same action I — I m . Using 
Eq. (13) we get a simple approximate relation for the 
difference between eigenvalues of neighboring modes [7] 



dm 



Hm+i — H m — u> m /k, 



(36) 



where uj m =u(I m ). 

The m-th eigenfunction (p m (z) between its turning 
points can be represented as [16] 



(37) 



where 



±i(kS m (z)+<f> b /2)) 



(38) 



S m {z) = / dzp m (z), p m {z) = ^2 [H m - U(z)}, 



and 



Qr, 



2-Kp m (z) ' 



(39) 
(40) 



The normal mode representation of the wave field in 
a range-independent waveguide is given by an expansion 
into a sum of eigenfunctions of the Sturm-Liouville eigen- 
value problem [8, 16] 

l^ + k 2 (H m -U) lfm = (32) 

with appropriate boundary conditions. The latter can 
be determined by reflection coefficients at the surface, 
V s = e % ^ s , and at the bottom, V b = e^ b . We assume 
that reflections at boundaries occur without energy loss 
and therefore <p s and 4>b are real constants. In particular, 
for a pressure release surface and rigid bottom we have 



B. Mode coupling equation 



In each cross-section of a waveguide the field can be 
decomposed into a sum of the local modes, that is the 
modes of the reference waveguide (see Sec. II C 2) cor- 
responding to this particular range. Since at different 
range points we deal with, generally, different reference 
waveguides, the eigenfunctions and eigenvalues become 
functions of range. The mode representation of the wave 
field in a range-dependent waveguide has the form 

u(r,z) = ^a m {r)tp m (r,z), (41) 

m 



■5 



where we have emphasized the r dependence of the eigen- 
function. By substituting this expression in Eq. (3) and 
using the orthogonality condition (33) we get the mode 
coupling equation 



da Tt 



dr 



/d 
dz —^-<p m . (42) 



For a point source defined by Eq. (10) it should be solved 
with an initial condition 



a m (0) = tp m (0,z o ). 



(43) 



For later convenience we present the derivative dip m /dr 
as an expansion 



dr 



(44) 



Substituting this in Eq. (42) yields 
da„ 



dr 



(45) 

v^O 



It can be shown that 



H„ 



1 f , dU 

—- / dz ip m ^-¥rn+v (46) 

- Hm+v J or 



This formula is derived by differentiating Eq. (32) for 
the reference waveguide with respect to r and exploiting 
the normalization condition (33). A detailed derivation 
of a more general result for the wave field governed by 
the Hclmholtz equation see in Ref. [8]. 



1. Matrix elements 

In the WKB approximation the matrix element of a 
smooth function F(z) [16] 



/ 



dz ip m (z)F(z)ip„(z) 



(47) 



with m,n » 1 can be expressed through coefficients of 
expansion (18). Indeed, take F(z) whose characteristic 
scale is much greater than the wavelength 2n/k. Then 
the matrix elements are very small unless |m — v\ << to. 
Making use of Eqs. (37)-(39) rewrite Eq. (47) in the 
form 

F m , m +Am = 2 J dz F{z)Q 2 m {z) 



X COS [k {S m -\-Am{z 

) - S m (z))] . 
From Eqs. (36) and (39) it follows that 



k (S m+A m(z) - S m {z)) = uj m Am 

J z„ 



dz 

m Pm(z)' 



(48) 



(49) 



According to Eq. (7) dz/p m = dr and the integral on the 
right can be considered as a shift along the r-axis between 
a minimum and a current point of a ray trajectory with 
action I m . In the range- independent waveguide u m r = 
9 and we get a desired connection between the matrix 
elements and the coefficients of Eq. (18) 



m,m+Am 



-!- / d9 F(z(I m ,9)) cos {Am 9) 
27T Jo 



-F\Am\(Im)- 



(50) 



According to this result Eq. (46) in the high frequency 
approximation reduces to 



B,, 



k 



2voj, 



-V M (I m ,r). 



(51) 



C. Adiabatic approximation for mode amplitudes 



In the adiabatic approximation the right hand side of 
the mode coupling equation (45) is assumed to be negli- 
gible [7, 8] and we get 



t (r) =a m {V)e- lkS « H ™ {r ' )dr ' ■ 



(52) 



The first order correction can be obtained using a simple 
perturbation theory. Consider the case when only one 
mode with number m is excited at r = 0, that is 



a m (0) = 5 V 



(53) 



Substituting Eq. (52) with a m (0) defined by Eq. (53) 
into the right hand side of Eq. (45) we find that for 

to 7^ mo 



(54) 



where 



'0 



dr'B m , mo _ m (r')e 4 ( m - m °)^'^o(-")^"_ (55) 



At high frequencies, when Eq. ( 51) is valid, q m can be 
presented in the form 



2(to — too) 



/' 

Jo 



^ T y, ,( r ' J )J(m-m ) fa' u mQ (r")dr" 



The condition 



\q m \ « f 



(56) 
(57) 
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(or its analog for the Hclmholtz equation) is traditionally 
considered as a starting point for studying the applica- 
bility of adiabatic approximation [7, 8]. In particular it 
turns out that the inequality (57) requires that 



D/L « 1, 



(58) 



where L is the characteristic scale of the horizontal range 
dependence. Another criterion involving the dependence 
on a frequency of the propagating wave has the form 
[8, 17] 



kD 2 /L « 1. 



(59) 



There are less general (but less restrictive and therefore 
more useful) criteria that can be obtained for waveguides 
with some special properties [8]. 



IV. PROJECTION OF THE RAY 
REPRESENTATION ONTO NORMAL MODES 

A. Mode amplitude as a function of ray 
parameters 

In this section we show that the mode amplitude in 
a range-dependent waveguide can be approximately ex- 
pressed through parameters of ray trajectories. In more 
details this issue has been considered in Rcfs. [1, 2]. Our 
starting point is an equation 



a m (r) = J dz v(r, z)(p m (r, z) 



(60) 



that follows from the mode orthogonality condition (33). 
Substituting u(r, z) from Eq. (8) and using the WKB 
approximation for <p m (r, z) we get a m = a+ + a~ with 



<„ = J dz B{r,z,a)e ik ^ r ' z '°\ 



where a denotes plus or minus, 

B = e -''/W2+«o-0t,/2 AQ m 

and 

<£ = S + aS m . 



(61) 



(62) 



(63) 



In the high frequency approximation (large k) integral 
(61) can be evaluated using a standard stationary phase 
technique [13]. This yields 



2tt 



£> exp 



k\d 2 <S>/dz 2 



ik$ + i^sgn (d 2 $/dz 



(64) 



with z st being a stationary phase point where the deriva- 
tive d$/dz vanishes. The function sgn(x) gives the sign 



of its argument. In classical mechanics it is well-known 
that the eikonal S considered as a function of range r, 
starting coordinate zo , and a final coordinate z obeys the 
relation [5] 



as 

dz 



= P, 



(65) 



where p is a momentum at range r. Combining this with 
Eq. (39) we get 



— = p + ap, n 



(66) 



This derivative vanishes at a point z such that a ray 
arriving at this point has an action 



and 



-sgn(p). 



(67) 



(68) 



The condition (67) singles out rays whose actions at the 
given range r are equal to the action of the m-th mode. 
In quantum theory a similar result was obtained in Ref. 
[4]- 

For the second derivative of $ we have an expression 



dz 2 



dp 

dz 

1 dH 

p dz 



_dU/dz 



uj_dl_ 

p dz ' 



(69) 



where H and / are considered as functions of the ray 
arrival depth at the given range r. From Eq. (64) it 
follows that 



cf n =A\kdI/dz\ 



-1/2 



exp (i/c$ — z7r/4CTSgn(9//i9z))| z=Zst . 



(70) 



To complete the calculation of the mode amplitude we 
should express the derivative dl /dz through parameters 
of rays. In the next subsection this is done for two im- 
portant examples of the starting field. 



B. Point source and source exciting a single mode 

In the case of a point source forming the starting field 
(10) we deal with a congruence of rays leaving point 
(0, zo). Labelling each ray by its starting momentum p 
we present the derivative dl/dz as (dl / dpo) / (dz / dpo) . 
Then insertion of Eq. (11) into Eq. (70) yields 

1 



y/2ir\dI/dpo\ 
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x exp (k (S + aS m ) - ifj.n/2 + iacj> b /2 + i((3 - 1)tt/4) , 

(71) 

with 



P = -<7Sgn(dI/dpo)sgn(dz/dpo). 



(72) 



Formula (71) accounts for a contribution from a ray sat- 
isfying the condition (67) to the ra-th mode. The mode 
amplitude is evaluated by summing up contributions 
from all such rays. Note that in the range-independent 
waveguide at any range point and for any mode there are 
exactly two rays satisfying Eq. (67) [1-3, 18]. They es- 
cape the source at grazing angles (equal in absolute value 
and opposite in sign) that coincides with grazing angles 



of the quasi-plane waves <p^ 



-ikH m r 



and <p m e 



ikH„ 



at 

depth z . This statement remains valid in an adiabatic 
waveguide as well. Notice that if rays are adiabatic, i.e. 
the term A in the Hamiltonian (21) is negligible, then 
Eq. (71) reduces to Eq. (52) with a m (0) defined by Eq. 
(43). A corresponding transformation is rather simple 
but somewhat lengthy and therefore we do not present it 
here. 

A similar result can be derived for a source exciting a 
single mode with m = m [1]. In this case we have two 
congruences of rays associated with functions <p mo . Ini- 
tial actions of all rays are equal to I mo and their starting 
depths are located between turning points of the mode. 
There are two rays escaping each point z within this 
interval of depths. These rays have starting momenta 



Pa = ±^2(H mo (0)-U(z )) 
and complex amplitudes [1, 19] 

A = Q m „(zq) p ± t ( Smo (z)+4> b /2) 

V\dz/dz \ 



(73) 



(74) 



Here we label rays by their starting depths. Substitut- 
ing Eq. (74) into Eq. (70) and representing 81/ dz as 
(dl / dzo) / (dz / dzo) yields 

Q mo (zv) 



^k\di/dz \ 



x exp (ik {S + aS m + f3 S mo (z ))) 



x exp (—ifj,n/2 + icr<t> b /2 + i/3tt/4) , 



(75) 



with 



= -asgn(dl/dz )sgn(dz/dz ), (3 = sgn(p ). (76) 

An important cautionary remark should be made. 
When using the stationary phase technique we assume 
that the integration in Eq. (61) goes over an interval 
exceeding 



Sz 



2tt 



h\d 2 ®\ 
K \~3z T \ 



(77) 



Therefore our approach requires that 

Sz « Az, 



(78) 



where Az is a difference between mode turning points. 
Denote by AI a spread of the ray action at the range of 
observation. A corresponding spread of the "energy" H 
in accord with Eqs. (13) and (24) is AH = 2-kAI/D. 
Assuming that the spread of ray depths at the range of 
observation has the same order of magnitude as Az and 
using Eq. (69) we obtain an order-of-magnitude estimate 



8 2 <S> 

dz 2 



1 AH 

p Az 



12vr AI 
p~~D~A~z' 



(79) 



where p is an rms value of momentum at range r. As it 
follows from Eq. (67) an initially excited mode breaks 
up into a group of approximately 



Am = kAI 

modes. Using a rough estimate pD/2 
ing Eqs. (77), (79), and (80) we get 



Sz: 



Az 



^/Am/2' 



(80) 

Az and combin- 



(81) 



The condition (78) is met only if Am >> 1, that is if 
many modes are excited. Thus, Eq. (75) is valid only in 
a strongly nonadiabatic environment and at long enough 
ranges where an initially excited mode breaks up into 
many modes. In contrast, a point source usually excites 
a large number of modes and therefore the condition (78) 
practically does not restrict applicability of Eq. (71). 



VARIATIONS OF ACTION VARIABLE AND 
MODE COUPLING 

Number of interacting modes and applicability 
of adiabatic approximation 



First, consider the case when only one mode is excited 
at r = 0. In spite of the remark made at the end of Sec. 
IV B the condition (67) even at short ranges properly in- 
dicates rays contributing to the given mode. Therefore 
Eq. (80) provides an estimate for the number of inter- 
acting modes valid at any range. Correspondingly, the 
applicability condition of adiabatic approximation can be 
formulated in the form 



AI « l/k. 



(82) 



The quantity AI on the left represents a difference be- 
tween initial and final actions along a typical ray. Sub- 
stituting Eq. (31) into the first of ray equations (23) 
yields 

AI = J(r) - I mo 
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1 OC „ r 



dr V v (I mo ,r)cos(u0) . 



(83) 



In the integrands on the right we neglect the deviation of 
action from its starting value. In the same approximation 
the second of ray equations (23) gives 



O + / u>(I mo ,r)dr. 



(84) 



Substituting this into Eq. (83) we see that Eq. (82) 
agrees with the condition (57). 

Since we consider a linear problem and the wave field 
can be always represented as a superposition of normal 
modes the criterion (82) remains valid for an arbitrary 
source. But in the case of a point source we have an ex- 
plicit expression for the mode amplitude (Eq. (71)) which 
can be used at both short ranges where the adiabatic ap- 
proximation is still valid and at long ranges where this 
approximation fails. Proceeding from this result we shall 
derive a more accurate criterion. It is natural to expect 
that the deviation from the adiabaticity first reveals it- 
self in a deviation of the phase k§ defined by Eqs. (63), 
(67), and (68) from its value obtained under assumption 
that A can be dropped. 

In order to estimate this phase deviation take a ray (we 
shall call it nonadiabatic) contributing to the m-th mode 
and compare its phase with that of a similar ray (adia- 
batic) whose trajectory satisfies Eqs. (23) with A = 0. 
Both rays escape the same point source and at range r 
have the same value of action I ~ I m . Formally, we shall 
assume that A is proportional to some small parameter e. 
Our purpose is to evaluate the difference between phases 
of the rays up to terms of order 0(e). We shall use the 
symbols 6S, 61, 69, and 6z to denote the difference in 
eikonals, action variables, angle variables, and vertical 
coordinates of the nonadiabatic and adiabatic rays, re- 
spectively. 

Compare eikonals (see Eq. (25)) of our rays. From 
Eqs. (15) it follows that 

G(zq, I m + SI Q ,0) - G(z ,I m ,0) = 9(0)61(0), (85) 



and 

G(z + 5z(r),I m ,r) - G(z,I m ,r) = p m Sz(r). 
Making use of ray equations (23) yields 



(86) 



6 f 1 (Id6 - Hdr) 

J rn 



ri ' dH 

I m dS0 + Side - —SIdr 
ol 



= I m (59(r) - 60(0)) . 
Combining Eqs. (86) and (87) we find 



6S = p m 6z( 



r)- f 



Adr. 



(87) 



(88) 



The integration goes along the adiabatic ray. We assume 
that the nonadiabaticity is so week that the constants 
/i, (3, and a for both rays are the same. According to 
Eq. (39) the difference between S m corresponding to our 
rays can be estimated as 5S m = p m 5z(r). Taking into 
account Eq. (68) we finally arrive at 



8<f> 



Adr. 



(89) 



A slightly different derivation of Eq. (89) is given in 
Rcf. [3]. Thus, applicability of adiabatic approximation 
requires that 



k I Adr 



< 7T. 



(90) 



Using an explicit expression for A given by Eq. (31) and 
approximating 9 by expression presented in Eq. (84) it is 
not difficult to show that the condition (90) agrees with 
Eqs. (56) and (57). 

The criterion (82) is much more convenient for practi- 
cal applications than (57) and (90). Indeed, evaluation 
of AI can be performed using a standard ray code with- 
out exploiting formula (83). To find action / at the given 
range r one should (i) compute the ray parameters p and 
z at this range and (ii) using them as initial conditions 
evaluate (with the same ray code) integral (12) over the 
ray cycle in a corresponding reference waveguide. 



B. Examples 

To verify and illustrate the above results we have com- 
puted wave fields in two range-dependent hydroacoustic 
waveguides. This is done using the code MMPE [20] 
originally created for solving the wide angle parabolic 
equation. It has been slightly modified to use the stan- 
dard parabolic equation approximation. All numerical 
results presented in this section have been obtained for 
monochromatic wave fields at a carrier frequency of 200 
Hz. 

In our first example the sound speed field is taken in 
the form c(r, z) = c(z) + 6c(r, z). A range independent 
constituent 



c(z) = c (l + e (e 2 ^-^' B - 2(z - z a )/B - 1 

' ' (91) 

with Co = 1.5 km/s, B = 1 km, z a = —1 km, and 
e = 0.0057 represents the so-called Munk profile widely 
used to study sound propagation in deep sea [7, 11]. 
We consider a strong range-dependent perturbation mod- 
elling a synoptic eddy 

Mr , 2) = C2ex p(-<!lZlH£-(i^# N ), (92) 



where 



Ar 2 



Az(r) = Az c — Az v exp 



Az(r) 2 

(r - r v ) 2 
Ar v 



(93) 



9 



The following values of parameters have been selected: 
c 2 = -0.01 km/s, r 2 = 300 km, z 2 = -1 km, Ar = 80 
km, Az c = 0.5 km, Az v = 0.25 km, r v = 320 km, Ar v = 
20 km. The isolines of the total sound field are shown in 
Fig. 1. 



1.5078 -~ ^ 




200 400 600 

r, km 



FIG. 1: Sound speed field modelling a synoptic eddy in a deep 
sea. Contour labels are of sound speed in km/s. 

The upper panel in Fig. 2 shows the deviation of 
mode amplitude at range r = 600 km from its start- 
ing value at r = 0. The deviation is taken relative to 
the rms mode amplitude at r = 0. It is seen that only 
modes with m < 15, whose amplitudes remain practi- 
cally unchanged, are adiabatic. This fact agrees with 
predictions following from Eqs. (82) and (90). The mid- 
dle and lower panels in Fig. 2 present values of 5$ and 
kAI = k(I(r) — I m ) computed along ray paths escap- 
ing the point source with starting values of action vari- 
ables equal to I m (to = 0, ... , 100). For each mode there 
are two such rays with launch angles equal in absolute 
value and opposite in sign (see comment after Eq. (68)). 
Solid (dashed) curves in both panels correspond to rays 
starting upward (downward). Consistent with our ex- 
pectation, although 8<fr/ (2tt) and k(I — I m ) arc not close 
they have the same order of magnitude and both are very 
small for adiabatic modes. 

The upper panel in Fig. 3 shows the range depen- 
dences of mode amplitudes in a situation when a single 
mode with m = 60 is excited at r = 0. Due to scatter- 
ing at the eddy the 60-th mode breaks up into a group 
of modes. In order to predict the width of this group 
using Eq. (80) we have traced a fan of 40 rays with ini- 
tial actions equal to Iqo- The rays start from 20 points 
uniformly sampling the depth interval between the turn- 
ing points of mode 60. There arc two rays with start- 
ing momenta defined by Eq. (73) escaping each point. 
The lower panel of Fig. 3 presents the deviation of ac- 
tion from its starting value as a function of range. It is 
clearly seen that in agreement with Eq. (80) kAI with 
AI = max \ I(r) — /go I representing the spread of actions 
at the given range, properly predicts the width of the 
group of excited modes. 

In another example of the range-dependent waveguide 




10 20 30 40 50 60 70 80 00 10O 

mode number, m 



FIG. 2: Modes excited by a point source at a carrier frequency 
200 Hz. The source is set at z = — 1 km. Upper panel: Nor- 
malized differences between mode amplitudes at ranges r — 
and r = 600 km; Aa m — a m (600 km) — a m (0), a rms is a rms 
mode amplitude at r — 0. Middle panel: phase variations i5$ 
due to nonadiabaticity predicted by Eq. (89) at r = 600 km 
for rays starting upward (solid line) and downward (dashed 
line) . Lower panel: deviations of ray actions from their start- 
ing values Im at r = 600 km computed for rays starting up- 
ward (solid line) and downward (dashed line). 




100 200 300 400 500 




FIG. 3: Upper panel: range dependences of mode amplitudes 
(shown by gradations of gray) for the case when only one 
mode with m = 60 is excited at r = 0. Lower panel: deviation 
of the action variable from its starting value computed for a 
fan of 40 rays starting from different depth points with actions 
1(0) = /«>. 

the sound speed field 

c(r, z) = ^pci(z) + ^c 2 (z), (94) 

with R = 200 km being the waveguide length, represents 
a linear superposition of two profiles C\(z) and c 2 (z). An 
evolution of the sound speed profile with range is shown 
in Fig. 4. In this waveguide the mode coupling is rather 
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strong and values of <5<I> and kAI (not shown) at range 
200 km are on the order of unity for all modes. To check 
the applicability of estimate (80) we, once again, have 
considered the case when only the 60-th mode is excited 
at r = 0. Figure 5 again demonstrates that kAI provides 
a good estimate for the width of a group of normal modes 
into which the original mode breaks up. A new feature 
of function AI(r) = max|7(r) — Iqo\ absent in Fig. 3 is 
an appearance of narrow spots in both panels at ranges 
of about 47, 94, 141, and 188 km. From the viewpoint of 
the ray-based approach this phenomenon is related to the 
fact that the sound speed field (94) is a linear function 
of range. Since U(r,z) ~ (c(r, z) — c )/c is an almost 
linear function of r, the derivative dU/dr is practically 
range-independent and the same is true of the coefficients 
V v present in Eq. (83). At not very long ranges 9 k 9 a + 
Lu(Ieo, 0)r and \I — Iqo\ in this approximation vanishes at 
ranges equal to integer multiples of the cycle length of 
the 60-th mode. The latter is about 47 km. 




sound speed (km/s) 



FIG. 4: Sound speed profiles at range points equally spaced 
within an interval from r — to r — 200 km are plotted 
(from left to right) with a sound speed offset of 0.4 km/s. 
The leftmost and rightmost curves represent functions ci(z) 
and ci[z) in Eq. (94), respectively. 
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VI. CONCLUSION 



In this paper we have considered a ray-based method 
representing an analog of the geometrical optics for 
modes. In the scope of this approach the mode ampli- 
tude is expressed through solutions of the ray equations. 
It turns out that the m-th mode is formed by contribu- 
tions from rays whose action variables at the range of 
observation, up to a multiplicative constant, are equal to 
m. Since the ray method remains valid in the nonadia- 
batic waveguide, the relatively simple formulas connect- 
ing rays and modes provide a convenient tool for study- 
ing the applicability of adiabatic approximation. In the 
present paper we have used this option. 

Our simplest criterion is given by Eq. (82) that im- 
poses limitations on the range variation of the action vari- 
able. This equation establishes connection between the 
validity of adiabatic approximations for rays and modes. 
It is interesting that when the condition (82) fails, the 
quantity kAI estimates a number of interacting modes. 
For the case of a point source we have derived a more ac- 
curate criterion (90) whose implementation, however, re- 
quires more detailed calculations. Nevertheless, it should 
be emphasized that (i) the left hand sides of both Eqs. 
(82) and (90) can be evaluated using a standard ray trac- 
ing and (ii) both criteria agree with the condition (57) 
derived from the mode coupling equation. Note also that 
unlike conditions (58) and (59), our criteria account for 
accumulation of errors with range. 

Since the left hand sides of Eqs. (82) and (90) are pro- 
portional to the carrier frequency it is clear that the lower 
the frequency, the better these inequalities are satisfied 
and, hence, the wider the validity region of the adiabatic 
approximation. The same conclusion follows from Eq. 
(59). 

Finally, notice that all our results have been obtained 
in the high frequency approximation which, naturally, 
restricts their generality. In contrast, the use of the 
parabolic equation approximation is not principal. The 
results can be easily generalized to the case when the 
wave field is governed by the Helmholtz equation. Al- 
though we have considered sound waves propagating 
in underwater acoustic waveguides the formulas derived 
here can be applied for description of wave propagation 
in different waveguide media and for analysis of a quan- 
tum particle oscillating in a potential well with time- 
dependent parameters. 

This work was supported by the Russian Foundation 
for Basic Research under Grant No. 03-02-17246. 



FIG. 5: The same as in Fig. 3 but constructed for a waveguide 
with the sound speed field defined by Eq. (94) and shown in 
Fig. 4. 
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